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ABSTRACT 

The energy spectra of galactic cosmic rays carry fundamental information regarding 
their origin and propagation. These spectra, when measured near Earth, are signifi- 
cantly affected by the solar magnetic field. A comprehensive description of the cosmic 
radiation must therefore include the transport and modulation of cosmic rays inside 
the heliosphere. During the end of the last decade the Sun underwent a peculiarly long 
quiet phase well suited to study modulation processes. In this paper we present proton 
spectra measured from July 2006 to December 2009 by PAMELA. The large collected 
statistics of protons allowed the time variation to be followed on a nearly monthly basis 
down to 400 MV. Data are compared with a state-of-the-art three-dimensional model 
of solar modulation. 

Subject headings: Cosmic rays; Heliosphere; Solar cycle; Solar modulation 

1. Introduction 

Protons are the most abundant species in cosmic rays, representing about 90% of the total 
flux. They are believed to be mainly of primordial origin, accelerated during supernovae explosions 
and confined within our Galaxy, at least up to PeV energies. During propagation through the 
interstellar medium cosmic rays undergo interactions with gas atoms and loose energy significantly 
modifying their injection spectra and composition. 

Beside the effect of propagation in the Galaxy, cosmic-ray particles reaching the Earth are 
affected by the solar wind and the solar magnetic field. A continuous flow of plasma emanates 
from the solar corona and extends out far beyond Pluto, up to the heliopause where the interstellar 
medium is encountered. The solar wind travels at supersonic speeds of about 350-650 km s^ 1 up to 
the termination shock where it slows down to about 100-150 km s _1 to form the inner heliosheath. 
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The plasma carries the solar magnetic field out into the heliosphere to form the heliospheric 
magnetic field (HMF). When interstellar cosmic rays enter the heliosphere, they interact with the 
solar wind and the HMF that modify their energy spectra. These modifications are directly related 
to solar activity and the total effect is called solar modulation. Solar activity varies strongly with 
time, rising from a minimum level when the Sun is quiet (with cosmic rays then having a maximum 
intensity at Earth) to a maximum period (when cosmic rays have a minimum intensity) and then 
returning to a new minimum to repeat the cycle. This solar cycle spans approximately 11 years. 
Solar activity modulates cosmic rays with the same cycle as a function of position in the heliosphere. 
Protons with rigidities up to at least 30 GV are affected and the effect becomes progressively larger 
as the rigidity decreases (e.g. iPotgieterl (|201ll )). At each solar activity maximum the polarity 



of the solar magnetic field reverses. When the solar magnetic field lines point outward in the 
Northern hemisphere and inward in the Southern hemisphere, the Sun is said to be in an A>0 cycle, 
whereas during an A<0 cycle the solar magnetic field lines point inward (outward) in the Northern 
(Southern) hemisphere. Direct cosmic-ray mea surements above a few MeV have been pe r forme d 
from the inner to the outer heliosphere (see e.g. iHeber and Potgieterl (|2006l ); IWebber et al.l (|2012i )) 
and particularly at the Earth where solar modulation effects are at their largest. 

Understanding the origin and propagation of cosmic rays requires knowledge of the cosmic-ray 
energy spectra in the interstellar medium, i.e. uninfluenced by the Sun's magnet ic field. These 



spectr a, especially for protons, can be inferred from gamma-ray data (e.g. see iNeronov et al 



(|2012j )). However, u ncertainties are large and direct cosmic-ray measurements are preferable. Voy- 
ager space missions ([Decker et al.ll2005l : lRichardson et al.ll2008l ) have passed the termination shock 
and may soon start sampling the cosmic radiation in the local interstellar medium. They will 
provide fundamental data to understand the basic features of the heliosphere. However, these 
cosmic-ray measurements are limited to the very low energy part of the spectra (from tens to a 
few 100 Me V), with galactic p rotons, helium and oxygen strongly contaminated by the anomalous 
component ( Stone et al.ll2008l ). Higher energies are needed to fully address fundamental questions 
in cosmic-ray physics. Moreover, the energy spectra of rare components like antiparticles may 
yield indirect information regarding new physics such as dark matter. For example, t he antiproton 



flux b e tween few hundred MeV- few GeV can be used (e.g. ICerdeno et al.l (|2012j ); lAsano et al 
(|2012j ); iGarnv et al.l (|2012j )) to constrain models with light (~ 10 GeV) dark matter candidates put 



forward to interpret recent measurement s by direct-detection experiments (jBernabei et al. 



Aalseth et al. 



2011 



2000 



Angloher et al.l 120121 ) . However, solar modulation has a relevant effect on the 



cosmic-ray spectra at these energies that needs to be disentangled to allow a comprehensive picture 
to emerge. 

Precise measurements of cosmic-ray spectra over a wide rigidity range, from a few hundred 
MV to tens of GV over an extended period of time (covering a significant fraction of the solar cycle) 
can be used to study the effect of solar modulation in greater detail, including the convective and 
adiabatic cooling effect of the expanding solar wind and the diffusive and particle drift effects of the 
turbulent HMF. The most recent solar minimum activity period (2006-2009) and the consequent 
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minimum modulation conditions for cosmic rays were unusual. It was expected that the new activity 
cycle would begin early in 2008. Instead minimum modulation conditions continued until the end 
of 2009 when the highest levels o f cosm ic rays were recorded at Earth since the beginning of the 
space age (see e.g. iMewaldt et al.l (|2010l )). This period of prolonged solar minimum activity is well 
suited to study the various modulation processes that affect the propagation of galactic cosmic rays 
inside the heliosphere. 

PAMELA is a satellite-borne instrument designed for cosmic-ray antimatter studies. The 
instrument is flying on-board the Russian Resurs-DKl satellite since June 2006, in a semi-pola r 
near-Earth or bit. PAMELA has pro vided important results on the antiproton (jAdriani et al. I l2009ah 
and positron (jAdriani et al.ll2009bl ) galactic abundances. The high-resol ution magnetic spect rom- 
eter allowed hydrogen and helium spectral measurements up to 1.2 TV (jAdriani et al.ll201ll ). the 



highest limit ever achieved by this kind of experiment. We present here results on the long-term 
variations in the absolute flux of hydrogen nuclei, measured down to 400 MV. 

The analysis, described in section[3l is based on data collected between July 2006 and December 
2009. The analyzed period covers the unusually long most recent solar minimum. The large proton 
statistics collected by the instrument allowed the proton flux to be measured for each Carrington 
rotation^] (hereafter the flux measured on this time basis will be referred as Carrington flux). 



The instrument 



The PAMELA instrument, shown schematically in Figure [U comprises a number of highly 
redundant detectors, capable of identifying particles by providing charge (Ze), mass, rigidity (R = 
pc/Ze, with p momentum of the particle and c speed of light) and velocity (/3) measurements over 
a very wide energy range. Multiple detectors are placed around a permanent magnet with a silicon 
micro-strip tracking system, providing charge and track deflection (d = ±1/R) information. A 
scintillator system (TOF: SI, S2 and S3) provides the trigger for event acquisition, the time-of- 
flight measurement and a further charge information. A silicon-tungsten calorimeter is used to 
perform hadron/lepton separation; a scintillator (S4) and a neutron detector placed at the bottom 
of the apparatus increase the rejection power at high energy. An anti-coincidence system (AC: 
CAS, CAT and CARD) registers hits produced by particles enterin g from outside the in strument 
acceptance. A detailed description of the apparatus can be found in lPicozza et al.1 (120071 ) . 



1 Mean synodic rotational period of the Sun surface, corresponding to about 27.28 days ( Carrington 18631 ) 



- 5 - 



TOF (S1) 



spectrometer 

tracking 
system 
(6 planes} 



geometric 
acceptance 




proton 



antiproton 



Fig. 1. — Schematic view of the PAMELA apparatus. 



3. Analysis 

3.1. Event selection 

In order to give a valid trigger for data acquisition, a particle must cross the S2 and S3 
scintillators, located above and below the tracking system, respectively. The requirement of a hit 
on S3 sets a limit on the minimum detectable rigidity, which for vertical protons is on average 
400 MV, i.e. ~ 80 MeV of kinetic energy. 

Clean events are selected by requiring a precision track in the spectrometer, a trajectory con- 
tained within the instrument acceptance and no hits in the CARD and CAT scintillators. Protons 
are identified among all the selected particles by requiring a positive track curvature and ionization 
energy losses (dE/dX) in the tracking system consistent with a hydrogen nucleus. The dE/dX selec- 
tion included deuterons, which in this work have been treated as being protons. Galactic particles 
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have been selected by requiring that the measured rigidity was a factor 1.3 above the geomagnetic 
vertical cutoff, evaluated using the Stormer approximation, at a given orbital position. The satellite 
orbit has an inclination of 70°, which allows galactic particles to be studied down to the minimum 
detectable rigidity. 

T he analysis procedur e was similar to previous work on the high energy proton and helium 
fluxes ( Adriani et al.l I2OI1I ). The major differences concern track-quality requirements and the 



fiducial volume definition. The procedure used in the proton and helium study aimed to maximise 
the spectrometer performance and minimize systematic uncertainties in the flux normalization, in 
order to perform a precise spectral measurement up to the highest achievable energy. The study 
presented in this paper aims to trace time variations in the low energy particle flux, with the finest 
possible time binning. Hence, high-performing tracking was less important than high statistics. 
For this reason the tracking requirements were relaxed (at least 3 hits, instead of 4, on both X and 
Y view, a track lever-arm of at least 4 silicon planes in the tracker) and the fiducial volume was 
enlarged (bounded 1.5 mm from the magnet cavity walls, instead of 7 mm), allowing to collect more 
particles with only a small degradation in the measurement precision in the energy range under 
study. 



3.2. Flux normalization 



The absolute proton flux in the spectrometer was obtained by dividing the measured energy 
spectrum by the acquisition time, the geometrical acceptance and the selection efficiency 



3.2.1. Acquisition time 

Flux evaluation was performed as a function of time, binning the data according to the Car- 
rington rotations. The time period studied corresponds to Carrington rotation numbers 2045-2092. 
The acquisition time was evaluated for each energy bin as the total live time spent above the geo- 
magnetic cutoff, outside the South Atlantic Anomaly and during periods when no solar events had 
occurred (hence December 2006, when a large solar event took place, was excluded from the data). 
This work is based on data collected between July 2006 and December 2009. The total live time 
was about 5xl0 7 s (~ 590 days) above ~ 20 GV, reducing to about 13% of this value at 0.4 GV, 
due to the relatively short time spent by the satellite at low geomagnetic locations. The accuracy 
of the live time determination was cross-checked by comparing different clocks available in flight, 
which showed a relative difference of less than 0.2 %. 
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3.2.2. Geometrical factor 

The geometrical factor for this analysis (G =19.9 cm 2 sr) is determined by the requirement of 
containment within the fiducial volume: a reconstructed track had to be 1.5 mm away from the 
magnet walls and the TOF-scintillator edges. The geometrical factor is constant above 1 GV and 
slowly decreases by ~ 2% to lower energies, where the curvature of the proton particle trajectory 
is no longer negligible. 

The value of G has been evaluated with two independent numerical methods, which provided 
values consistent within 0.1%. 



3.2.3. Selection efficiency 

The PAMELA instrument has been designed to provide redundant information on event topol- 
ogy and particle characteristics. This feature allowed the efficiencies of each detector to be measured 
directly from flight data. Due to the large proton statistics, detector behavior could be monitored 
over very short time scales. For this analysis, efficiencies were evaluated independently for each 
time interval, in order to follow any possible time variation of detector responses. 

The major time dependent effect was the sudden failure of some front-end chips in the tracking 
system. Chip failures started to occur few months after the launch of the satellite, initially with 
an almost constant rate. In the following years the rate of failures significantly decreased and 
the tracking system reached an almost stable condition near the end of 2009. The net effect 
was a progressive reduction of the tracking efficiency, since the number of hits available for track 
reconstruction decreased. No degradation in the signal-to-noise ratio and spatial resolution was 
observed. 



The approach described in lAdriani et al.l (|201ll ) for the evaluation of the tracking efficiency 



was used: a high energy average value was experimentally obtained using non-interacting minimum 
ionizing particles in the calorimeter from flight data and an energy dependent correction factor was 
evaluated using Monte Carlo simulation. The average tracking efficiency was measured using a test 
sample of non-interacting relativistic protons in the imaging calorimeter. The resulting efficiency 
varied almost monotonically from ~ 94% during the first months to ~ 20% in the final months. 
The efficiency variation was correctly reproduced by the simulation, which included a map of dead 
channels as a function of time. Since the requirement of non-interacting minimum ionizing particles 
in the calorimeter selected essentially protons with rigidities > 2 GV, the determination of the 
tracking efficiency at lower rigidities was performed using the simulation. The simulation results 
were normalized to experimental values at high rigidities after performing an identical selection. 
This was required because, due to the relatively poor tracking resolution of the calorimeter, the 
average efficiency values obtained using the calorimeter selection were different from the average 
efficiency of an isotropically-distributed particle sample (as was the case for triggered events). This 
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normalization was applied to the simulation estimated efficiency that accounted for the full energy 
dependence of the tracking efficiency and was used for the flux estimation. When comparing the 
average experimental and simulated efficiencies, differences of less than 4% were found, except 
during the last few months 2009 when the tracking efficiency fell below 40% and, consequently, 
the discrepancy became more significant, increasing up to 13%. These differences were included 
in the systematic uncertainties o f the measuremen ts. As a further check, two different sim ulation 
packages were used: GEANT3 ferun et al.1 1 19941 ) and GEANT4 (jAgostinelli et al.1 12003I ). The 
resulting efficiencies differed by ~ 1.5% and this was included in the systematic uncertainty. 

TOF and AC selection efficiencies were evaluated relative to the tracking selection. For well 
contained events, the TOF and AC responses were not correlated and the selection efficiencies could 
be measured independently. At low rigidities TOF and AC efficiency measurements were affected by 
background contamination in the test samples of low energy secondary particles produced by high 
energy primaries interacting with the apparatus. This contamination resulted in an underestimation 
of the efficiency, since background interactions were often associated with multiple hits in both the 
TOF and AC scintillators. An upper limit on this contamination was determined by excluding the 
TOF and AC conditions from the test sample selection, thereby adding background events into the 
test sample. The difference between the results was found to be ~ 6% in the lowest-energy bin. 
This value was conservatively treated as an asymmetric systematic uncertainty on the absolute 
flux. 



3.3. Spectral unfolding and other corrections 

The proton spectrum measured in the spectrometer was corrected to the top of the payload by 
accounting for the effects of rigidity displacement due to finite spectrometer resolution and particle 
slow-down by ionization energy loss, the latter being more relevant at low energies. The correction 
was applied by m eans of an unfolding procedure, following the Bayesian approach proposed by 



D'Agostinil (j 19951 ). The detector response matrix was obtained from the simulation. The resulting 
spectrum was corrected for the contamination of locally produced secondary particles and for 
attenuation due to particles interacting with the satellite and the instrument itself. In order to 
estimate the residual contamination of locally produced secondary particles, simulations of protons 
and helium nuclei isotropically hitting the PAMELA payload within a 2n solid angle were carried 
out. The estimated background was found to be less than 6% at 0.4 GV and inversely proportional 
to energy, becoming negligible above 1 GV. Flux attenuation, which was slightly energy dependent 
and amounted to 5-6%, accounted for both particle loss by non-elastic interactions and loss and 
gain of particles from within the acceptance due to elastic scattering. 

Finally, in order to correct for systematic time-dependent variations of the Carrington fluxes 
(see next section), these were normalized, at hi gh-energy (30- 5 GeV ), to the proton flux measured 



with the more stringent tracking conditions of lAdriani et al.l (|201ll ) and averaged over the period 
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Fig. 2. — Carrington proton fluxes shown relative to the flux obtained as in lAdriani et al.l ( 20 111 ) 
(value averaged over the period July 2006 - March 2008) for (open circles) stringent tracking criteria 
and (filled circles) relaxed tracking criteria. Two energy intervals are presented. 



July 2006 - March 20ofl 



3.4. Systematic uncertainties 

By relaxing the tracking selection the proton statistics could be significantly increased, allowing 
the proton flux to be evaluated on a Carrington rotation basis up to the end of 2009. This approach 
also led to an increase of the systematic uncertainties on the measured absolute fluxes. 

The reliability of the flux determination was verified by studying the time variation of the flux 
at high energy (above 30 GeV) since no significant modulation effects are expected in this energy 
region. As a reference, we considered the fluxes obtained with the stringent tracking conditions 
and averaged over the period July 2006 - March 2008. Figure [2] shows the time dependence of the 
ratio between the Carrington fluxes and the reference flux for two high-energy intervals at about 
30 GV and 47 GV. Two different selection criteria were considered: stringent (open circles, as in 



An improved use of the simu lation resulted in sl ightly higher selection efficiencies than previously estimated. 
Therefore, the proton spectrum of lAdriani et al.l q201ir ) should be lowered by 3.2%. Furthermore, the overall system- 
atic uncertainty is slightly larger at low energies (4.5% at 1 GV). 
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Adriani et al.l (|201ll )) and relaxed tracking selections (filled circles, this work). 

Comparing the upper and lower panels of Figure [2] a correlation between the time profiles at 
different energies can be noticed, which indicates the pres ence of non-negligibl e systematic effects. 
The Carrington fluxes based on the analysis described in lAdriani et al.l ( 2011 ) fluctuate around a 
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Fig. 3. — The time profile of the proton flux for three energy intervals (from top to bottom): 
0.597-0.625 GeV, 0.99-1.04 GeV and 2.99-3.13 GeV. The error bars are statistical and the shaded 
area represents the quadratic sum of all systematic uncertainties. 

mean value consistent with the quoted ~ 4% systematic uncertainty, while the fluxes for this analysis 
show a decreasing trend, with a relative variation over the whole time interval of ~10%. This time- 
dependent systematic effect was corrected as discussed in the previous section. The assumption 
behind this correction implies that PAMELA results are not sensitive to solar modulation induced 
flux variations smaller than about 5%. 



Further contributions to the total uncertainty resulted from the unfolding procedure (1-5%) 
and from the correction to the top of the payload (!%)• 



4. Results 

Figure [3] shows the time profile of the proton flux for three energy intervals; the error bars 
are statistical, while the shaded area represents the quadratic sum of all systematic uncertainties. 
The systematic uncertainty is dominated by the tracking efficiency and reaches a minimum above 
a few tens of GeV and for the first Carrington rotations, when the tracking efficiency is maximum. 
Additional time and energy-dependent contributions to the overall uncertainty have been introduced 
in order to account for both the reduced accuracy of the simulation in reproducing the spectrometer 
behavior as the detection efficiency decreases in time, and the larger uncertainty on the simulated 
efficiency at low energy. 

Figure [4] shows the evolution of the Carrington proton energy spectra from July 2006 (blue 
curve) to December 2009 (red curve). Around 100 MeV proton intensities increased by ~ 20% over 
twelve months, from December 2006 to December 2007, with a slightly larger increase of ~ 25% 
from December 2007 to December 2008. For 2009, when solar minimum conditions continued to 
prevail throughout the heliosphere, intensities increased by almost 40%. From these data it is clear 
that protons (evident across the greater part of the energy spectrum) reached maximum intensities 
at the end of 2009. 



5. Data Interpretation and Discussion 

Figure [Hand Table [T] present the galactic proton spectra measured by PAMELA over four time 
periods. These data illustrate how the proton spectra evolved from late 2006 to late 2009. The 
latter represents the highest proton flux observed by PAMELA. In Figure [5] the PAMELA proton 
spectra are overlaid with the corresponding computed spectra (solid lines). A full three-dimensional 
model was used to compute the differential intensity of cosmic-ray protons at Earth from 1 MeV 



to 30 GeV. The model is based on the numerical solution of the Parker transport equation ([Parker 



19651 ). including all four major modulation mechanisms: convection, particle drift caused by gra- 
dients and curvatures in the HMF, diffusion described by a full 3D diffusion tensor and adiabatic 
energy changes. For these calculations a local interstellar spectrum ( LIS) was assumed (shown a s 



a solid black line in Figure [5]). The proton LIS was based on that by lLangner &; Potgieterl (J2004) 
modified at high energies to match PAMELA observations. It follows from the comparison between 
the model and observations that the proton spectra became progressively softer, as expected, from 
2006 but requiring larger diffusion coefficients at lower energies than anticipated. At the same 
time full particle drifts occurred but with the modulation drift effects somewhat obscured by the 
effective diffusion. This caused the effect of the HMF wavy current sheet to be less well correlated 
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Fig. 4. — The evolution of the proton energy spectrum as particle intensities approached the period 
of minimum solar activity, from July 2006 (violet), to December 2009 (red). The region between 
the blue and red curves indicates the spread in proton fluxes during this time. 



to increases in cosmic-ray intensities at Earth. Details concerning the modulation model and theo- 
retica l assumptions and implications, will be published in an accompanying paper (jPotgieter et al 
20121 ). 
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Table 1. Hydrogen flux measured by PAMELA over five time periods. The first and second errors represent the statistical 

and systematic uncertainties, respectively. 



Kinetic Energy Flux 

(GeV) (m 2 s sr McV)- 1 

2006/11/13 - 2006/12/4 2007/11/30-2007/12/27 2008/11/19-2008/12/15 2009/12/6-2010/1/1 






082 


- 0.095 


(0 


721 


± 





032 


+0.061-. 
-0.100.' 


(1 


045 


± 





031 


+0.090 
-0.102 


(1 


191 


± 





035 


+0.105-, 
-0.104/ 


(1 


730 


± 





143 


+0.266-, 
-0.163-* 





095 


- 0.10 


(0 


759 


± 





036 


+0.063 \ 
-0.096/ 


(1 


111 


± 





030 


+0.094 
-0.098 


(1 


332 


± 





034 


+0.116-, 

-0.100/ 


(1 


722 


± 





156 


+0.263-, 
-0.160-* 





10 - 


0.11 


(0 


818 


± 





030 


+0.067-, 
-0.093/ 


(1 


151 


± 





027 


+0.097-, 
-0.094,1 


(1 


401 


± 





031 


+0.121-, 
-0.097,* 


(2 


016 


± 





134 


+0.307-, 
-0.158,* 





11 - 


0.12 


(o 


893 


± 





034 


+0.072-, 
-0.089/ 


(1 


230 


± 





027 


+0.102-, 
-0. 091.> 


(1 


477 


± 





031 


+0.126-, 
-0.094/ 


(2 


210 


± 





151 


+0.336-, 
-0.157,* 





12 - 


0.13 


(o 


928 


± 





030 


+0.074-. 
-0.086,* 


(1 


334 


± 





026 


+0.110-, 
-0.088,> 


(1 


543 


± 





029 


+0.131 \ 
-0.091/ 


(2 


162 


± 





129 


+0.327-, 
-0.155,* 





13 - 


0.15 


(o 


960 


± 





028 


+0.076-, 
-0.084.> 


(1 


411 


± 





025 


+0.115-, 
-0.086> 


(1 


607 


± 





027 


+0.135 -, 
-0.089/ 


(2 


256 


± 





123 


+0.341 -, 

-0.154,* 





15 - 


0.16 


(1 


030 


± 





028 


+0.081 -, 
-0.083 -* 


(1 


453 


± 





023 


+0.117-, 
-0.085.* 


(1 


675 


± 





026 


+0.139-, 
-0.087/ 


(2 


447 


± 





122 


+0.368-, 
-0.153/ 





16 - 


0.17 


(1 


085 


± 





027 


+0.084-, 
-0.081 - 1 


(1 


477 


± 





022 


+0.118-, 
-0.083.* 


(1 


727 


± 





025 


+0.142 -, 
-0.085 > 


(2 


359 


± 





112 


+0.354-, 
-0.152,* 





17 - 


0.19 


(1 


123 


± 





025 


+0.086-, 
-0.079.> 


(1 


545 


± 





021 


+0.122-, 
-0.081.* 


(1 


805 


± 





024 


+0.147-, 
-0.084,1 


(2 


277 


± 





102 


+0.341-, 
-0.151,* 





19 - 


0.21 


(1 


164 


± 





024 


+0.088-, 
-0.078.> 


(1 


613 


± 





020 


+0.126-, 
-0.080.* 


(1 


865 


± 





023 


+0.150-, 
-0.082,' 


(2 


423 


± 





097 


+0.362 \ 
-0.150,* 





21 - 


0.22 


(1 


194 


± 





022 


+0.090-, 
-0.076.> 


(1 


634 


± 





019 


+0.126-, 
-0.078,1 


(1 


899 


± 





021 


+0.152-, 
-0.081 -* 
+0.156-, 
-0.080 .* 


(2 


592 


± 





093 


+0.386-, 
-0.149,* 





22 - 


0.24 


(1 


222 


± 





022 


+0.091 -, 
-0.075,* 


(1 


673 


± 





018 


+0.128-, 
-0.077 ,* 


(1 


968 


± 





021 


(2 


647 


± 





091 


+0.393-, 
-0.149,1 





24 - 


0.26 


(1 


278 


± 





021 


+0.094-, 
-0.074,1 


(1 


732 


± 





018 


+0.131-, 
-0.076,1 


(2 


031 


± 





020 


+0.159-, 
-0.079 .* 


(2 


682 


± 





088 


+0.397-, 
-0.148,* 





26 - 


0.29 


(1 


359 


± 





020 


+0.098-, 
-0.073-1 


(1 


775 


± 





017 


+0.133-, 
-0.075 1 


(2 


047 


± 





019 


+0.159-, 
-0.078/ 
+0.155 -, 
-0.077/ 


(2 


684 


± 





082 


+0.396-, 
-0.148,* 





29 - 


0.31 


(1 


371 


± 





020 


+0.098-, 
-0.072,1 


(1 


788 


± 





016 


+0.132-, 
-0.074.* 


(2 


025 


± 





018 


(2 


603 


± 





077 


+0.383-, 
-0.147/ 





31 - 


0.34 


(1 


380 


± 





019 


+0.098-, 
-0.071 I 


(1 


818 


± 





016 


+0.133-, 
-0.073J 


(2 


040 


± 





017 


+0.155 -, 
-0.076 ) 


(2 


656 


± 





075 


+0.390-, 
-0.147,* 





34 - 


0.36 


(1 


381 


± 





018 


+0.097-, 


(1 


822 


± 





015 


+0.132-, 
-0.072 ,* 


(2 


036 


± 





017 


+0.153-, 
-0.075.* 


(2 


633 


± 





072 


+0.385-, 
-0.146,* 





36 - 


0.39 


(1 


398 


± 





017 


~0 097 
-0.069/ 


(1 


816 


± 





014 


+0. 130 \ 
-0.072 ) 


(2 


019 


± 





016 


+ 0. 150 -, 
-0.074,* 


(2 


504 


± 





067 


+0.366 \ 
-0.146,* 





39 - 


0.43 


(1 


410 


± 





016 


+0.096-, 
-0.068/ 


(1 


778 


± 





013 


+0.126-, 
-0.071.* 


(1 


983 


± 





015 


+0.146-, 
-0.074.* 


(2 


415 


± 





061 


+0.351-, 
-0.146,* 





43 - 


0.46 


(1 


396 


± 





016 


+0.094-, 
-0.068/ 


(1 


754 


± 





013 


+0.123-, 
-0.070.* 


(1 


959 


± 





014 


+0.143-, 
-0.073/ 


(2 


495 


± 





060 


+0.362-, 
-0.145.* 





46 - 


0.50 


(1 


391 


± 





015 


+0.093-, 
-0.067,> 


(1 


715 


± 





012 


+0.119-, 
-0.069.* 


(1 


941 


± 





014 


+0.140-, 
-0.072,* 


(2 


383 


± 





056 


+0.345-, 
-0.145 J 





50 - 


0.54 


(1 


349 


± 





014 


+0.089-, 
-0.066,1 


(1 


677 


± 





012 


+0.115-, 
-0.068.* 


(1 


844 


± 





013 


+0.132-, 
-0.071 .* 


(2 


200 


± 





052 


+0.318-, 
-0.144,* 





54 - 


0.58 


(1 


328 


± 





014 


+0.086-, 
-0.065/ 


(1 


643 


± 





011 


+0.111-1 

-0.068.* 
+0.105-, 
-0.067.) 


(1 


789 


± 





012 


+0.126-, 
-0.071 > 


(2 


213 


± 





050 


+0.319-, 
-0.144.* 





58 - 


0.62 


(1 


300 


± 





009 


+0.084-, 
-0.064,1 


(1 


577 


± 





010 


(1 


731 


± 





013 


+0.121-, 
-0.070 -* 


(2 


107 


± 





032 


+0.303-, 
-0.144/ 





62 - 


0.67 


(1 


254 


± 





008 


+0.080-, 
-0.064,1 


(1 


527 


± 





009 


+0.101-, 
-0.066.* 


(1 


641 


± 





012 


+0.114-, 
-0.069 .* 


(1 


966 


± 





030 


+0.282-, 
-0.143-* 





67 - 


0.72 


(1 


220 


± 





008 


+0.077-, 
-0.063,> 


(1 


468 


± 





009 


+0.096-, 
-0.065.) 


(1 


587 


± 





011 


+0.109-, 
-0.068.* 


(1 


870 


± 





028 


+0.268-, 
-0.143,* 





72 - 


0.77 


(1 


173 


± 





007 


+0.073-, 
-0.062,1 


(1 


403 


± 





008 


+0.091-, 
-0.065.* 


(1 


510 


± 





011 


+0.102-, 
-0.068 -> 


(1 


769 


± 





026 


+0.252-, 
-0.143-* 





77 - 


0.83 


(1 


125 


± 





007 


+0.069-, 
-0.061/ 


(1 


341 


± 





007 


+0.086-, 
-0.064.* 


(1 


431 


± 





010 


+0.096-, 
-0.067-* 


(1 


710 


± 





025 


+0.243-, 
-0.142,* 





83 - 


0.89 


(1 


064 


± 





006 


+0.064-, 
-0.060 J 


(1 


272 


± 





007 


+0.080-, 
-0.063.1 


(1 


357 


± 





009 


+0.090-, 
-0.066 .* 


(1 


550 


± 





022 


+0.220-, 
-0.142,1 





89 - 


0.96 


(1 


003 


± 





006 


+0.060-, 
-0.060 J 


(1 


197 


± 





007 


+0.075-, 
-0.062.* 


(1 


273 


± 





009 


+0.084-, 
-0.066 .* 


(1 


442 


± 





021 


+0.204-, 
-0.142,* 





96 - 


1.02 


(o 


964 


± 





005 


+0.057-, 
-0.059,> 


(1 


135 


± 





006 


+0.070-, 
-0.062.* 


(1 


192 


± 





008 


+0.077-, 
-0.065-* 


(1 


382 


± 





020 


+0.195-, 
-0.141,* 


1 


02 - 


1.09 


(o 


919 


± 





005 


+0.053-, 
-0.058,> 


(1 


066 


± 





006 


+0.065-, 
-0.061.* 


(1 


127 


± 





007 


+0.072-, 
-0.064,* 


(1 


294 


± 





018 


+0.183-, 
-0.141,* 


1 


09 - 


1.17 


(o 


860 


± 





005 


+0.049-, 
-0.057,' 


(o 


993 


± 





005 


+0.060-, 
-0.060 ) 


(1 


041 


± 





007 


+0.066-, 
-0.064.1 


(1 


178 


± 





017 


+0.166-, 
-0.141.* 
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Table 1 — Continued 



Kinetic Energy 




Flux 




(GcV) 




(m 2 s sr MeV)" 1 




2006/11/13 - 2006/12/4 


2007/11/30-2007/12/27 


2008/11/19-2008/12/15 


2009/12/6-2010/1/1 



1.17 - 


1.25 


(0.813 


± 





004 


± 


0.046) 




(0.934 


± 





005 + 056*1 








(0.983 


± 





006 -1- 062*1 








(1.105 


± 





016 + 155*1 








1.25 - 


1.34 


(0.753 


± 





004 


± 


0.042) 




(0.860 


± 





005 ± 0.051) 








(0.903 


± 





006 ± 0.056) 








(1.030 


± 





014 ± 0.144) 








1.34 - 


1.42 


(0.708 


± 


o 


004 


± 


0.039) 




(0.813 


± 


o 


004 + 047*1 








(0.845 


± 


o 


006 + 052*1 








(0.931 


± 


o 


01 3 + 1 30*1 








1.42 - 


1.52 


(0.650 


± 


o 


004 


± 


0.035) 




(0.748 


± 


o 


004 + 043*1 








(0.782 


± 


o 


005 -1- 048*1 








(0.862 


± 


o 


012 + 120*1 








1.52 - 


1.62 


(0.607 


± 


o 


003 


± 


0.033) 




(0.693 


± 


o 


004 + 039*1 








(0.719 


± 


o 


005 -1- 044*1 








(0.789 


± 


o 


Oil + 1 1 0*1 








1.62 - 


1.72 


(0.562 


± 


o 


003 


± 


0.030) 




(0.640 


± 


o 


003 + 036*1 








(0.670 


± 


o 


004 + 040*1 








(0.720 


± 


o 


010 + 100*1 








1.72 - 


1.83 


(0.520 


± 


o 


003 


± 


0.027) 




(0.590 


± 


o 


003 ± 0.033) 








(0.612 


± 


o 


004 ± 0.036) 








(0.675 


± 


o 


010 ± 0.094) 








1.83 - 


1.95 


(0.481 


± 


o 


003 


± 


0.025) 




(0.540 


± 


o 


003 + 030*1 








(0.562 


± 


o 


004 -1- 033*1 








(0.605 


± 


o 


009 + 084*1 

v *y j •.' i_ y ) . v/v " i ) 








1.95 - 


2.07 


(0.441 


± 


o 


002 


± 


0.023) 




(0.501 


± 


o 


003 + 027*1 








(0.512 


± 


o 


004 + 030*1 








(0.547 


± 


o 


008 + 076*1 

\J\JyJ 1_ \ r . \ > t \j j 








2.07 - 


2.20 


(0.412 


± 


o 


002 


± 


0.021) 




(0.450 


± 


o 


002 ± 0.024) 








(0.461 


± 


o 


003 ± 0.027) 








(0.499 


± 


o 


008 ± 0.069) 








2.20 - 


2.33 


(0.374 


± 


o 


002 


± 


0.019) 




(0.409 


± 


o 


002 ± 0.022) 








(0.428 


± 


o 


003 ± 0.024) 








(0.450 


± 


o 


007 ± 0.062) 








2.33 - 


2.48 


(0.347 


± 


o 


002 


± 


0.017) 




(0.377 


± 


o 


002 + 020*1 








(0.388 


± 


o 


003 + 022*1 








(0.411 


± 


o 


006 + 057*1 

v *y j \i i_ y i . v / : ) \ i 








2.48 - 


2.62 


(0.315 


± 


o 


002 


± 


0.015) 




(0.348 


± 


o 


002 ± 0.018) 








(0.356 


± 


o 


003 ± 0.020) 








(0.390 


± 


o 


006 ± 0.054) 








2.62 - 


2.78 


(0.286 


± 


o 


002 


± 


0.014) 




(0.318 


± 


o 


002 ± 0.016) 








(0.322 


± 


o 


002 ± 0.018) 








(0.334 


± 


o 


005 ± 0.046) 








2.78 - 


2.94 


(0.262 


± 


o 


002 


± 


0.012) 




(0.288 


± 


o 


002 ± 0.015) 








(0.291 


± 


o 


002 ± 0.016) 








(0.303 


± 


o 


005 ± 0.042) 








2.94 - 


3.12 


(0.239 


± 


o 


001 


± 


0.011) 




(0.258 


± 


o 


002 ± 0.013) 








(0.262 


± 


o 


002 ± 0.014) 








(0.276 


± 


o 


005 ± 0.038) 








3.12 - 


3.30 


(0.215 


± 


o 


001 


± 


0.010) 




(0.234 


± 


o 


001 ± 0.012) 








(0.239 


± 


o 


002 ± 0.013) 








(0.251 


± 


o 


004 ± 0.034) 








3.30 - 


3.49 


(0.195 


± 





001 


± 


0.009) 




(0.210 


± 





001 ± 0.010) 








(0.215 


± 





002 ± 0.012) 








(0.223 


± 





004 ± 0.031) 








3.49 - 


3.69 


(0.175 


± 





001 


± 


0.008) 




(0.189 


± 





001 ± 0.009) 








(0.194 


± 





002 ± 0.010) 








(0.204 


± 





004 ± 0.028) 








3.69 - 


4.12 


(151.4 


± 





7 ± 


6.9) x 10" 


3 


(163.3 


± 





7565 ± 8.050) 


x 


10 


-3 


(164.7 


± 





9853 ± 8.790) 


X 


10 


-3 


(0.172 


± 





002 ± 0.023) 








4.12 - 


4.59 


(124.4 


± 





6 ± 


5.6) x 10" 


3 


(132.0 


± 





6391 ± 6.446) 


x 


10 


-3 


(134.3 


± 





8364 ± 7.113) 


X 


10 


-3 


(0.141 


± 





002 ± 0.019) 








4.59 - 


5.11 


(99.15 


± 





5 ± 


4.4) x 10" 


3 


(106.8 


± 





5402 ± 5.168) 


x 


10 


-3 


(109.1 


± 





7083 ± 5.731) 


X 


10 


-3 


(0.113 


± 





002 ± 0.015) 








5.11 - 


5.68 


(81.38 


± 





4 ± 


3.6) x 10" 


3 


(85.34 


± 





4537 ± 4.093) 


X 


10 


-3 


(86.97 


± 





5944 ± 4.536) 


X 


10 


-3 


(0.091 


± 





001 ± 0.012) 








5.68 - 


6.30 


(65.27 


± 





4 ± 


2.9) x 10" 


■3 


(68.61 


± 





3820 ± 3.262) 


X 


10 


-3 


(69.55 


± 





4994 ± 3.601) 


X 


10 


-3 


(0.069 


± 





001 ± 0.009) 








6.30 - 


6.99 


(52.12 


± 





3 ± 


2.3) x 10" 


-3 


(54.86 


± 





3211 ± 2.587) 


x 


10 


-3 


(55.60 


± 





4197 ± 2.859) 


X 


10 


-3 


(56.69 


± 





9533 ± 7.693) 


X 


io- 


-3 


6.99 - 


7.74 


(41.66 


± 





3 ± 


1. 


8) x 10- 


-3 


(43.65 


± 





2691 ± 2.042) 


x 


10 


-3 


(44.15 


± 





3514 ± 2.255) 


X 


10 


-3 


(44.91 


± 





7967 ± 6.088) 


X 


io- 


-3 


7.74 - 


8.57 


(33.18 


± 





2 ± 


1.4) x 10- 


-3 


(34.31 


± 





2241 ± 1.593) 


x 


10 


-3 


(34.65 


± 





2923 ± 1.759) 


X 


10 


-3 


(34.70 


± 





6572 ± 4.700) 


X 


io- 


-3 


8.57 - 


9.48 


(26.26 


± 





2 ± 


1.1) x 10- 


-3 


(26.80 


± 





1856 ± 1.235) 


X 


10 


-3 


(27.12 


± 





2424 ± 1.368) 


X 


10 


-3 


(27.44 


± 





5473 ± 3.713) 


X 


io- 


-3 


9.48 - 


10.48 


(208.8 


± 


1 


5 ± 


8.8) x 10- 


-4 


(214.3 


± 


1 


556 ± 9.872) x 10" 


4 


(21.32 


± 





2014 ± 1.075) 


X 


10 


-3 


(21.77 


± 





4564 ± 2.945) 


X 


io- 


-3 


10.48 


- 11.57 


(164.0 


± 


1 


2 ± 


6.9) x 10- 


4 


(167.1 


± 


1 


286 ± 7.697) x 10" 


4 


(167.7 


± 


1 


672 ± 8.456) x 10" 


4 


(16.86 


± 





3755 ± 2.281) 


X 


io- 


-3 


11.57 


- 12.77 


(128.9 


± 


1 


± 


5.5) x 10- 


4 


(130.2 


± 


1 


061 ± 5.998) x 10" 


4 


(134.6 


± 


1 


399 ± 6.785) x 10" 


4 


(13.03 


± 





3083 ± 1.763) 


X 


io- 


-3 


12.77 


- 14.09 


(100.1 


± 





8 ± 


4.2) x 10- 


4 


(103.4 


± 





8845 ± 4.762) 


x 


10 


-4 


(102.8 


± 


1 


143 ± 5.182) x 10" 


4 


(10.22 


± 





2554 ± 1.384) 


X 


io- 


-3 


14.09 


- 15.54 


(78.59 


± 





7 ± 


3.3) x 10- 


4 


(79.72 


± 





7250 ± 3.672) 


x 


10 


-4 


(78.55 


± 





9319 ± 3.960) 


X 


10 


-4 


(8.339 


± 





2150 ± 1.129) 


X 


io- 


-3 



Table 1 — Continued 



Kinetic Energy 




Flux 




(GcV) 




(m 2 s sr McV) -1 




2006/11/13 - 2006/12/4 


2007/11/30-2007/12/27 


2008/11/19-2008/12/15 


2009/12/6-2010/1/1 



15.54 - 


17.12 


(62.71 


± 





6 


± 


2.7) 


X 


io- 


-4 


(64.08 


± 


0.6064 ± 2.952) x 10~ 4 


(62.57 


± 


0.7736 ± 3.154) x IO" 4 


(62 


71 


± 


1.735 ± 8.486) x 


10~ 4 


17.12 - 


18.86 


(47.92 


± 





5 


± 


2.0) 


X 


io- 


-4 


(49.08 


± 


0.4937 ± 2.261) x 10~ 4 


(49.15 


± 


0.6357 ± 2.478) x IO" 4 


(47 


39 


± 


1.400 ± 6.413) x 


10~ 4 


18.86 - 


20.76 


(37.16 


± 





4 


± 


1.6) 


X 


io- 


-4 


(37.96 


± 


0.4066 ± 1.748) > 


< io- 4 


(38.06 


± 


0.5230 ± 1.919) x IO" 4 


(38 


42 


± 


1.181 ± 5.200) x 


io- 4 


20.76 - 


22.85 


(28.92 


± 





3 


± 


1.2) 


X 


io- 


-4 


(29.32 


± 


0.3398 ± 1.351) > 


( 10~ 4 


(28.47 


± 


0.4282 ± 1.435) x IO" 4 


(28 


76 


± 


0.9688 ± 3.892) > 


io- 4 


22.85 - 


25.15 


(225.8 


± 


2 


6 


± 


9.6) 


X 


io- 


-5 


(23.07 


± 


0.2879 ± 1.063) > 


( io- 4 


(22.61 


± 


0.3642 ± 1.140) x 10~ 4 


(21 


59 


± 


0.8009 ± 2.922) > 


io~ 4 


25.15 - 


27.66 


(171.3 


± 


2 


2 


± 


7.2) 


X 


io- 


-5 


(173.8 


± 


2.389 ± 8.006) x 


10" 5 


(173.5 


± 


3.049 ± 8.745) x IO" 5 


(18 


35 


± 


0.7055 ± 2.484) x IO" 4 


27.66 - 


30.42 


(138.0 


± 


1 


9 


± 


5.8) 


X 


io- 


-5 


(137.7 


± 


2.031 ± 6.344) x 


IO" 6 


(138.8 


± 


2.603 ± 6.997) x IO" 5 


(14 


03 


± 


0.5883 ± 1.898) x IO" 4 


30.42 - 


33.44 


(103.8 


± 


1 


5 


± 


4.4) 


X 


io- 


-5 


(105.6 


± 


1.699 ± 4.862) x 


IO" 5 


(101.5 


± 


2.127 ± 5.118) x IO" 5 


(10 


95 


± 


0.4965 ± 1.482) x 10~ 4 


33.44 - 


36.75 


(81.42 


± 


1 


3 


± 


3.4) 


X 


io- 


-5 


(80.24 


± 


1.416 ± 3.696) x 


10~ s 


(81.51 


± 


1.820 ± 4.109) x IO" 5 


(76 


20 


± 


3.955 ± 10.31) x 


IO" 5 


36.75 - 


40.39 


(63.41 


± 


1 


1 


± 


2.7) 


X 


io- 


-5 


(62.97 


± 


1.199 ± 2.900) x 


IO" 6 


(62.22 


± 


1.520 ± 3.137) x IO" 5 


(61 


51 


± 


3.396 ± 8.323) x 


10~ 5 


40.39 - 


44.37 


(47.29 


± 





9 


± 


2.0) 


X 


io- 




(49.30 


± 


1.014 ± 2.271) x 


IO" 8 


(49.86 


± 


1.300 ± 2.514) x 10~ 5 


(47 


38 


± 


2.847 ± 6.412) x 


io- 5 


44.37 - 


48.74 


(35.22 


± 





8 




1.5) 


X 


io- 




(33.67 




0.801 ± 1.551) x 


IO" 5 


(35.50 


± 


1.049 ± 1.790) x 10~ 5 


(36 


99 




2.404 ± 5.006) x 


IO" 5 



